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Abstract 

In this paper, we study the problem of moving n sensors on a line to form a barrier coverage 
of a specified segment of the line such that the maximum moving distance of the sensors is 
minimized. Previously, it was an open question whether this problem on sensors with arbitrary 
sensing ranges is solvable in polynomial time. We settle this open question positively by giving 
, an 0{n^ logn) time algorithm. For the special case when all sensors have the same-size sensing 

range, the previously best solution takes 0{n?) time. We present an O(nlogn) time algorithm 
for this case; further, if all sensors are initially located on the coverage segment, our algorithm 
takes 0{n) time. Also, we extend our techniques to the cycle version of the problem where the 
barrier coverage is for a simple cycle and the sensors are allowed to move only along the cycle. 
■ For sensors with the same-size sensing range, we solve the cycle version in 0{n) time, improving 

the previously best O(n^) time solution. 

o[ 

^ ■ 1 Introduction 

O ■ A Wireless Sensor Netvi^ork (WSN) uses a large number of sensors to monitor some surrounding 

environmental phenomena [1]. Each sensor is equipped with a sensing device with limited battery- 
^ . supplied energy. The sensors process data obtained and forward the data to a base station. Intrusion 

I detection and border surveillance constitute a major application category for WSNs. A main goal 

CD ■ of these applications is to detect intruders as they cross the boundary of a region or domain. 

For example, research efforts were made to extend the scalability of WSNs to the monitoring of 
. international borders \10\ lllj . Unlike the traditional full coverage \13\ \TE[ [T9] which requires an 

Q ■ entire target region to be covered by the sensors, the barrier coverage [21 [3l O O [11] only seeks to 

cover the perimeter of the region to ensure that any intruders are detected as they cross the region 
border. Since barrier coverage requires fewer sensors, it is often preferable to full coverage. Because 
sensors have limited battery-supplied energy, it is desired to minimize their movements. In this 
^ ■ paper, we study a one-dimensional barrier coverage problem where the barrier is for a (finite) line 

segment and the sensors are initially located on the line containing the barrier segment and allowed 
to move on the line. As discussed in the previous work [HJ |9l [16] and shown in this paper, barrier 
coverage even for 1-D domains poses some challenging algorithmic issues. Also, our 1-D solutions 
may be used in solving more general problems. For example, if the barrier is sought for a simple 
polygon, then we may consider each of its edges separately and apply our algorithms to each edge. 

In our problem, each sensor has a sensing range (or range for short) and we want to move the 
sensors to form a coverage for the barrier such that the maximum sensor movement is minimized. 
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We present efficient algoritlims for tliis problem, improving the previous work and settling an open 
question. Also, we extend our techniques to the cycle version where the barrier is for a simple cycle 
and the sensors are allowed to move only along the cycle. 

1.1 Problem Definitions, Previous Work, and Our Results 

Denote hy B = [0, L] the barrier that is a line segment from x = 0tox = L>0on the x-axis. A 
set S" = { si, S2, • • • , Sn} of n mobile sensors is initially located on the x-axis. Each sensor Si € S has 
a range > and is located at the coordinate Xi of the x-axis. We assume xi < X2 < • • • < Xn- If a 
sensor Si is at the position x' , then we say Si covers the interval [x' — ri,x' + r^], called the covering 
interval of Sj. Our problem is to find a set of destinations on the x-axis, {yi, 1/2, ■ ■ ■ , Un}, for the 
sensors (i.e., for each Si G S, move Sj from Xj to yi) such that each point on the barrier B is covered 
by at least one sensor and the maximum moving distance of the sensors (i.e., maxi<j<„{|xj — yi\}) 
is minimized. We call this problem the barrier coverage on a line segment, denoted by BCLS. We 
assume 2 • Y17=i — (otherwise, a barrier coverage for B is not possible). 

The decision version of BCLS is defined as follows. Given a value A > 0, determine whether 
there is a feasible solution in which the moving distance of each sensor is at most A. The decision 
version characterizes a problem model in which the sensors have a limited energy and we want to 
know whether their energy is sufficient to move and form a barrier coverage. 

If the ranges of all sensors are the same (i.e., the r^'s are all equal), then we call it the uniform 
case of BCLS. When the sensors have arbitrary ranges, we call it the general case. 

The BCLS problem has been studied before. The uniform case has been solved in O(n^) time [8]. 
An 0(n) time algorithm is also given in [8j for the decision version of the uniform case. However, 
it has been open whether the general case is solvable in polynomial time [8]. 

In this paper, we settle the open problem on the general BCLS, presenting an 0(n^ logn) time 
algorithm for it. We also solve the decision version of the general BCLS in O(nlogn) time. Since 
this is a basic problem on sensors and intervals and our algorithm is the first-known polynomial 
time solution for it, we expect our results and techniques to be useful for other related problems. 
Further, for the uniform case, we derive an 0{n log n) time algorithm, improving the previous O(n^) 
time solution [8]; if all sensors are initially on B, our algorithm runs in 0{n) time. 

In addition, we consider the simple cycle barrier coverage where the barrier is represented as 
a simple cycle and the n sensors are initially on the cycle and are allowed to move only along the 
cycle. The goal is to move the sensors to form a barrier coverage and minimize the maximum sensor 
movement. If all sensors have the same range, Mehrandish [15] gave an O(n^) time algorithm, and 
we present an 0(n) time solution in this paper. 

1.2 Related Work 

Besides the results mentioned above, an 0(n) time 2-approximation algorithm for the uniform 
BCLS was also given in [8J and a variation of the decision version of the general BCLS is shown to 
be NP-hard [8]. Additional results were also given in [8j for the case 2 • Y17=i^i ^ (although in 
this case, B cannot be entirely covered). 

Mehrandish et al. |151ll6j also considered the line segment barrier, but unlike the BCLS problem, 
they intended to use the minimum number of sensors to form a barrier coverage, which they proved 
to be NP-hard. But, if all sensors have the same range, polynomial time algorithms were possible 
[15\ [TU] . Another study of the line segment barrier [S] aimed to minimize the sum of the moving 
distances of all sensors; this problem is NP-hard [S], but is solvable in polynomial time when all 
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sensors have the same range [9]. In addition, Li et al. [12] considers the hnear coverage problem 
which aims to set an energy for each sensor to form a coverage such that the cost of all sensors is 
minimized. There [12], the sensors are not allowed to move, and the more energy a sensor has, the 
larger the covering range of the sensor and the larger the cost of the sensor. 

Bhattacharya et al. |2] studied a 2-D barrier coverage problem in which the barrier is a circle and 
the sensors, initially located inside the circle, are moved to the circle to form a coverage such that 
the maximum sensor movement is minimized; the ranges of the sensors are not explicitly specified 
but the destinations of the sensors are required to form a regular n-gon on the circle. Subsequent 
improvements of the results in [2] have been made [U [T7]. In addition, Bhattacharya et al. [2] 
presented some results on the corresponding min-sum problem version (minimizing the sum of the 
moving distances of all sensors); further improvement was also given in [H I17j. 

Some other barrier coverage problems have been studied. For example, Kumar et al. [lij pro- 
posed algorithms for determining whether a region is barrier covered after the sensors are deployed. 
They considered both the deterministic version (the sensors are deployed deterministically) and 
the randomized version (the sensors are deployed randomly), and aimed to determine a barrier 
coverage with high probability. Chen et al. [3J introduced a local barrier coverage problem in which 
individual sensors determine the barrier coverage locally. 

1.3 An Overview of Our Approaches 

For any problem we consider, let A* denote the maximum sensor movement in an optimal solution. 

For the uniform BCLS, as shown in [8], a useful property is that there always exists an order 
preserving optimal solution, i.e., the order of the sensors in the optimal solution is the same as that 
in the input. Based on this property, the previous O(n^) time algorithm [8] covers B from left to 
right; in each step, it picks the next sensor and re-balances the current maximum sensor movement. 
In this paper, we take a very different approach. With the order preserving property, we determine 
a set A of candidate values for A* with A* € A. Consequently, by using the decision algorithm, we 
can find A* in A. But, this approach may be inefficient since |A| = Q{n?). To reduce the running 
time, our strategy is not to compute the set A explicitly. Instead, we compute an element in A 
whenever we need it. A possible attempt would be to first find a sorted order for the elements of A 
or (implicitly) sort the elements of A, and then obtain A* by binary search. However, it seems not 
easy to (implicitly) sort the elements of A. Instead, based on several new observations, we manage 
to find a way to partition the elements of A into n sorted lists, each list containing 0(n) elements. 
Next, by using a technique called binary search on sorted arrays [5], we are able to find A* in A in 
0(n log n) time. For the special case when all sensors are initially located on B, a key observation 
we make is that A* is precisely the maximum value of the candidate set A. Although A = 0(n^), 
based on new observations, we show that its maximum value can be computed in 0{n) time. 

For the general BCLS, as indicated in [8], the order preserving property no longer holds. Con- 
sequently, our approach for the uniform case does not work. The main difficulty of this case is that 
we do not know the order of the sensors appeared in an optimal solution. Due to this difficulty, 
no polynomial time algorithm was known before for the general BCLS. To solve this problem, we 
first develop a greedy algorithm for the decision version of the general BCLS. After 0(n log n) time 
preprocessing, our decision algorithm takes 0{n) time for any value A. If A > A*, implying that 
there exists a feasible solution, then our decision algorithm can determine the order of sensors in 
a feasible solution for covering B. For the general BCLS, we seek to simulate the behavior of the 
decision algorithm on A = A*. Although we do not know the value A*, our algorithm determines 
the same sensor order as it would be obtained by the decision algorithm on the value A = A*. 
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To this end, each step of the algorithm uses our decision algorithm as a decision procedure. The 
idea is somewhat similar to parametric search [6l |T3], and here we "parameterize" our decision 
algorithm. However, we should point out a few differences. First, unlike the typical parametric 
search [6l[l4], our approach does not involve any parallel scheme and is practical. Second, normally, 
if a problem can be solved by parametric search, then there also exist other (simpler) polynomial 
time algorithms for the problem although they might be less efficient than the parametric search 
solution (e.g., the slope selection problem [7]). In contrast, for our general BCLS problem, so far 
we have not found any other (even straightforward) polynomial time algorithm. 

In addition, our 0{n) time algorithm for the simple cycle barrier coverage is a generalization of 
our approach for the special case of the uniform BCLS when all sensors are initially located on B. 

For ease of exposition, we assume that initially no two sensors are located at the same position 
(i.e., Xi 7^ Xj for any i ^ j), and the covering intervals of any two different sensors do not share a 
common endpoint. Our algorithms can be easily generalized to the general situation. 

The rest of the paper is organized as follows. In Section [21 we describe our algorithms for the 
general BCLS. In Section [21 we present our algorithms for the uniform BCLS. Our results for the 
simple cycle barrier coverage are discussed in Section [H 

2 The General Case of BCLS 

In this section, we present our algorithms for the general BCLS problem. Previously, it was an 
open problem whether the general BCLS can be solved in polynomial time. The main difficulty is 
that we do not know the order of the sensors in an optimal solution. Our main effort is for resolving 
this difficulty, and we derive an O(n^logn) time algorithm for the general BCLS. 

We first give our algorithm for the decision version (in Section 12. ip , which is crucial for solving 
the general BCLS (in Section 12. 2p that we refer to as the optimization version of the problem. 

For each sensor Si S S, we call the right (resp., left) endpoint of the covering interval of Si the 
right (resp., left) extension of Sj. Each of the right and left extensions of Si is an extension of Sj. 
Denote by p{x') the point on the x-axis whose coordinate is x' , and denote by p^(x') (resp., p~{x')) 
a point to the right (resp., left) of p(x') and infinitely close to p{x'). The concept p~^{x') and 
p~{x') is only used to explain the algorithms, and we never need to find such a point explicitly. 
Let A* denote the maximum sensor moving distance in an optimal solution for the optimization 
version of the general BCLS problem. Note that we can easily determine whether A* = 0, say, in 
0(n log n) time. Henceforth, we assume A* > 0. 

2.1 The Decision Version of the General BCLS 

Given any value A, the decision version is to determine whether there is a feasible solution in 
which the maximum sensor movement is at most A. Clearly, there is a feasible solution if and 
only if A > A*. We show that after 0(n log n) time preprocessing, for any A, we can determine 
whether A > A* in 0{n) time. We explore some properties of a feasible solution in Section [2.1.1^ 
describe our decision algorithm in Section 12.1.21 argue its correctness in Section 12.1.31 and discuss 
its implementation in Section 12.1.41 In Section 12.1.51 we show that by extending the algorithm, 
we can also determine whether A > A* in the same time bound; this is particularly useful to our 
optimization algorithm in Section [2. 2[ 
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2.1.1 Preliminaries 



By a sensor configuration, we refer to a specification of where each sensor G 5" is located. By this 
definition, the input is a configuration in which each sensor Sj is located at Xi. The displacement of 
a sensor in a configuration C is the distance between the position of the sensor in C and its original 
position in the input. A configuration C is a feasible solution for the distance A if the sensors in C 
form a barrier coverage of B (i.e., the union of the covering intervals of the sensors in C contains 
B) and the displacement of each sensor is at most A. In a feasible solution, a subset S" C 5 is 
called a solution set if the sensors in 5" form a barrier coverage; of course, S itself is also a solution 
set. A feasible solution may have multiple solution sets. A sensor Si in a solution set S' is said to 
be critical with respect to S' if Sj covers a point on B that is not covered by any other sensor in 
iS". If every sensor in S' is critical, then S' is called a critical set. 

Given any value A, if A > A*, then our decision algorithm will find a critical set and determine 
the order in which the sensors of the critical set will appear in a feasible solution for A. For the 
purpose of giving some intuition and later showing the correctness of our algorithm, we first explore 
some properties of a critical set. 

Consider a critical set S'^. For each sensor s G S"^, we call the set of points on B that are covered 
by s but not covered by any other sensor in S"^ the exclusive coverage of s. 

Observation 1 The exclusive coverage of each sensor in a critical set S'^ is a continuous portion 
of the barrier B. 

Proof: Assume to the contrary the exclusive coverage of a sensor s G 5"^ is not a continuous 
portion of B. Then there must be at least one sensor s' G 5"^ whose covering interval is between 
two consecutive continuous portions of the exclusive coverage of s. But that would mean s' is not 
critical since the covering interval of s' is contained in that of s. Hence, the observation holds. □ 
For a critical set S"^ in a feasible solution SOL, we define the cover order of the sensors in S"^ 
as the order of these sensors in SOL such that their exclusive coverages are from left to right. 

Observation 2 The cover order of the sensors of a critical set S'^ in a feasible solution SOL is 
consistent with the left-to-right order of the positions of these sensors in SOL. Further, the cover 
order is also consistent with the order of the right (resp., left) extensions of these sensors in SOL. 

Proof: Consider any two sensors Sj and Sj in S^ with ranges and rj, respectively. Without loss 
of generality, assume Sj is to the left of Sj in the cover order, i.e., the exclusive coverage of Sj is to 
the left of that of Sj in SOL. Let y,; and yj be the positions of Si and Sj in SOL, respectively. To 
prove the observation, it suffices to show yi < yj, yi + ri < yj + rj, and yi — ri < yj — rj. 

Let p he a point in the exclusive coverage of Sj . We also use p to denote its coordinate on the 
X-axis. Then p is not covered by Si, implying either p > -\- ri or p < yi — ri. But, the latter case 
cannot hold (otherwise, the exclusive coverage of Si would be to the right of that of Sj). Since p 
is covered by Sj, we have p < yj + rj. Therefore, yi + ri < p < yj + rj. By using a symmetric 
argument, we can also prove yi — ri < yj — rj (we omit the details). Clearly, the two inequalities 
Vi + fi < Vj + i^j and yi — ri < yj — rj imply yi < yj. The observation thus holds. □ 

An interval / of S is called a left-aligned interval if the left endpoint of J is at (i.e., I is of 
the form [0, a;'] or [0,.x')). A set of sensors is said to be in attached positions if the union of their 
covering intervals is a continuous interval of the x-axis whose length is equal to the sum of the 
lengths of these covering intervals. Two intervals of the x-axis overlap if they intersect each other 
(even at only one point). 
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2.1.2 The Algorithm Description 

Initially, we move all sensors of S to the right by the distance A, i.e., for each 1 < z < n, we move 
Si to the position x'^ = Xi + A. Let Cq denote the resulting configuration. Clearly, there is a feasible 
solution for A if and only if we can move the sensors in Cq to the left by at most 2A to form a 
coverage of B. Thus, henceforth we only need to consider moving the sensors to the left. Recall 
that we have assumed that the extensions of any two distinct sensors are different; hence in Cq, the 
extensions of all sensors are also different. 

Our algorithm takes a greedy approach. It seeks to find sensors to cover B from left to right, 
in at most n steps. If A > A*, the algorithm will end up with a critical set S'^ of sensors along with 
the destinations for all these sensors. In theory, the other sensors in \ S'^ can be anywhere such 
that their displacements are at most A; but in the solution found by our algorithm, they are at the 
same positions as in Cq. If a sensor is at the same position as in Cq, we say it stands still. 

In step i (initially, i = 1), using the configuration Cj_i produced in step i — 1 and based on 
certain criteria, we find a sensor Sg(j) and determine its destination yg(i), where g{i) is the index of 
the sensor in S and ^^(j) € [^^^(j) — 2A,x^j^j^]. We then move the sensor Sg(j) to ^^(j) to obtain a new 
configuration Cj from Cj_i (if ^^(j) = then we need not move Sg(j), and Cj is simply Cj_i). 

Let Ri = y g[i) + r g(^i-) (i.e., the right extension of Sg(j) in Ci). Assume Rq = 0. Let Si = Sj-i U {sg(j)} 
(5*0 = initially). We will show that the sensors in Si together cover the left-aligned interval [0, Ri\. 
If Ri > we have found a feasible solution with a critical set S^ = Si, and terminate the algorithm. 
Otherwise, we proceed to step i + 1. Further, it is possible that a desired sensor Sg(j) cannot be 
found, in which case we terminate the algorithm and report A < A*. Below we give the details, and 
in particular, discuss how to determine the sensor Sg(j) in each step. 

Before discussing the first step, we provide some intuition. Let Si consist of the sensors whose 
right extensions are at most in Cq. We claim that since L > 0, no sensor in Si can be in a critical 
set of a feasible solution if A* < A. Indeed, because all sensors have been moved to their rightmost 
possible positions in Cq, if no sensor in Si has a right extension at in Cq, then the claim trivially 
holds; otherwise, suppose sj is such a sensor. Assume to the contrary that sj is in a critical set S'^. 
Then |?(0) is the only point on B that can be covered by st- Since L > 0, there must be another 
sensor in S^ that also covers p(0) (otherwise, no sensor in S^ would cover the point p^(0)). Hence, 
St is not critical with respect to S'^, a contradiction. The claim thus follows. Therefore, we need 
not consider the sensors in Si since they do not help in forming a feasible solution. 

In step 1, we determine the sensor as follows. Define = {sj \ x'^ — rj < < x'j + rj} 

(see Fig. [1]), i.e., Su consists of all sensors covering the point p(0) in Cq except any sensor whose 
right extension is (but if the left extension of a sensor is 0, the sensor is included in Su). In other 
words, Su consists of all sensors covering the point ^"''(0) in Cq. If ^ 0, then we choose the 
sensor in Su whose right extension is the largest as Sg(i) (e.g., Si in Fig. [1]), and let yg(i) = ^^(i)- 
Note that since the extensions of all sensors in Cq are different, the sensor Sg(i) is unique. If = 0, 
then define as the set of sensors whose left extensions are larger than and at most 2 A (e.g., see 
Fig. [2]). If Si2 = 0, then we terminate the algorithm and report A < A*. Otherwise, we choose the 
sensor in whose right extension is the smallest as Sg(i) (e.g., Sj in Fig. [2|), and let yg(i) = ?^g(i) 
(i.e., the left extension of Sg(i) is at after it is moved to the destination yg(i))- 

If the algorithm is not terminated, then we move Sg(i) to yg[i), yielding a new configuration Ci. 
Let 5i = {sg(i)}, and Ri be the right extension of Sg(i) in Ci. If Ri > L, we have found a feasible 
solution Ci with the critical set ^i, and terminate the algorithm. Otherwise, we proceed to step 2. 

The general step is very similar to step 1. Consider step i for i > 1. We determine the sensor 
Sg(^i), as follows. Let Su be the set of sensors covering the point p+(i2j_i) in the configuration Ci-i. 
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Figure 1: The set 5*11 consists of the three sensors whose Figure 2: The set S12 consists of the three sensors whose 
covering intervals are shown, and Sg(i) is Sj. covering intervals are shown, and Sg(i) is s, if Sn = 0. 

If Sii 0, we choose the sensor in Sn with the largest right extension as Sg(j) and let y^^j) = x'^^^y 
Otherwise, let 5^2 be the set of sensors whose left extensions are larger than Ri-i and at most 
Ri-i + 2A. If Si2 = 0, we terminate the algorithm and report A < A*. Otherwise, we choose the 
sensor in Si2 with the smallest right extension as and let yg(j) = Ri-i + fg^iy If the algorithm 
is not terminated, we move s^j-j) to y^^j) and obtain a new configuration Cj. Let Si = Si-i U 
Let Ri be the right extension of s^^j) in Q. U Ri > L, we have found a feasible solution Q with 
the critical set Si and terminate the algorithm. Otherwise, we proceed to step i + 1. If the sensor 
Sp(j) is from Sn (resp., Si2), then we call it the type I (resp., type II) sensor. 

Since there are n sensors in S, the algorithm is terminated in at most n steps. This finishes the 
description of our algorithm. 

2.1.3 The Correctness of the Algorithm 

Based on the description of our algorithm, we have the following lemma. 

Lemma 1 At the end of step i, suppose the algorithm produces the set Si and the configuration 
Ci; then Si and Ci have the following properties. 

(a) Si consists of sensors that are type I or type II. 

(b) For each sensor Sgf^j^ G Si with 1 < j < i, if ■'^g{j) is of type I, then it stands still (i.e., its 

position in Ci is the same as that in Cq); otherwise, its left extension is at Rj-i, and SgQ^ 
and ^^(j-i) are in attached positions if j > 1. 

(c) The interval on B covered by the sensors in Si is [0, Ri] . 

(d) For each 1 < j <i, the right extension of Sg(^j) is larger than that of Sg^j^iy 

(e) For each 1 < j < i, Sg(j) is the only sensor in Si that covers the point p'^{Rj-i) (with Rq = 0). 

Proof: The first three properties are trivially true according to the algorithm description. 

For property (d), note that the right extension of Sg(j) (resp., SgQ_i-j) is Rj (resp., Rj-i). 
According to our algorithm, the sensor Sg^-^ covers the point p'^{Rj-i), implying that Rj > Rj-i- 
Hence, property (d) holds. 

For property (e), note that the sensor s^q) (which is determined in step j) always covers 
p'^{Rj-i). Consider any other sensor 5^(4) G Si. If t < j, then the right extension of s^^^) is at most 
Rj-i, and thus Sg(j) cannot cover p'^(Rj-i). If t > j, then depending on whether s^^j) G Sti or 
Sp(t) G St2, there are two cases. If G St2, then the left extension of Sg(^t) is Rt-i, which is larger 
than Rj-i, and thus cannot cover p'^{Rj^i) in Cj. Otherwise (i.e., Sg^^-^ G Sti), stands 
still. Assume to the contrary that 5^(4) covers p^{R.j^i) in Cj. Then Sg(^t) must have been in Sji in 
step J within the configuration Cj-i. This implies Sji / 0, s^^j^ G Sji, and s^^j^ stands still. Since 
Rt is the right extension of Sp(() and Rj is the right extension of by property (d), for t > j, 



we have Rt > Rj. Since Rt > Rj (i.e., the right extension of is smaller than that of Sg(^t)), the 
algorithm cannot choose from Sji in step j, which is a contradiction. Therefore, ■Sg(f) cannot 
cover the point Property (e) thus holds. □ 

At its termination, our algorithm either reports A > A* or A < A*. To argue the correctness 
of the algorithm, below we will show that if the algorithm reports A > A*, then indeed there is a 
feasible solution for A and our algorithm finds one; otherwise, there is no feasible solution for A. 

Suppose in step i, our algorithm reports A > A*. Then according to the algorithm, it must be 
Ri > L. By Lemma [UJ^c) and[T]J^e), Ci is a feasible solution and Si is a critical set. Further, by 
Lemma [T]||d) and Observation [2l the cover order of the sensors in Si is Sg(i), •5^(2), • • • , 

Next, we show that if the algorithm reports A < A*, then certainly there is no feasible solution 
for A. This is almost an immediate consequence of the following lemma. 

Lemma 2 Suppose S^ is the set of sensors in the configuration Ci whose right extensions are at 
most Ri. Then the interval [0,Ri] is the largest possible left-aligned interval that can be covered by 
the sensors of S[ such that the displacement of each sensor in S[ is at most A. 

Proof: In this proof, when we say an interval is covered by the sensors of S[, we mean (without 
explicitly stating) that the displacement of each sensor in S[ is at most A. 

We first prove a key claim: If C is a configuration for the sensors of S[ such that a left-aligned 
interval [0, x'] is covered by the sensors of S[, then there always exists a configuration C* for S[ in 
which the interval [0, x'] is still covered by the sensors of S[ and for each 1 < j < i, the position of 
the sensor Sp(j) in C* is yg(j), where g{j) and y^Q) are the values computed by our algorithm. 

As similar to our discussion in Section 12.1.11 the configuration C for S*,- always has a critical set 
for covering the interval [0, x']. Let Sc be such a critical set of C. 

We prove the claim by induction. We first show the base case: Suppose there is a configuration 
C for the sensors of S^ in which a left-aligned interval [0, x'] is covered by the sensors of S'^; then 
there is a configuration C( for the sensors of in which the interval [0, x'] is still covered by the 
sensors of S'^ and the position of the sensor in C[ is yg(i) ■ 

Let t = g{l). If the position of st in C is yt, then we are done (with C[ = C). Otherwise, let y[ 
be the position of st in C, with y[ ^ yt. Depending on st G or st G S12, there are two cases. 

• If St € then yt = x[. Since yt is the rightmost position to which the sensor st is allowed 
to move and y't^ yt, we have y't < yt- Depending on whether st is in the critical set Sc, there 
further are two subcases. 

If St ^ Sc, then by the definition of a critical set, the sensors in Sc form a coverage of [0, x'] 
regardless of where st is. If we move st to yt (and other sensors keep the same positions as in 
C) to obtain a new configuration C[, then the sensors of S'^ still form a coverage of [0, x']. 

If St G Sc, then because yt > y't, if we move st from y'^ to yt, st is moved to the right. Since 
st S when st is at yt, st still covers the point p(0). Thus, moving st from y't to yt does 
not cause sj to cover a smaller sub-interval of [0, x']. Hence, by moving st to yt, we obtain a 
new configuration C( in which the sensors of S'^ still form a coverage of [0,x']. 

• If St £ S12, then according to our algorithm, 5ii = in this case, and st is the sensor in 
with the smallest right extension in Cq. If st Sc, then by the same argument as above, we 
can obtain a configuration C( in which the interval [0, x'] is still covered by the sensors of S"^ 
and the position of the sensor st in C{ is yt. Below, we discuss the case when st € Sc- 
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Figure 3: Illustrating the switch operation on St and Sh'- The left figure is before the switch and the right one is 
after the switch. 

In Sc, some sensors must cover the point p{0) in C. Let S' be the set of sensors in Sc that 
cover p{0) in C. If st G S' , then it is easy to see that y[ < yt since yt is the rightmost position 
for St to cover p{0). In this case, again, by the same argument as above, we can always move 
St to the right from y'^ to yt to obtain a configuration C[ in which the interval [0, x'] is still 
covered by the sensors of S'^. Otherwise (i.e., st ^ S'), we show below that we can always 
move Sf to yt by switching the relative positions of St and some other sensors in Sq- 

An easy observation is that each sensor in S' must be in Su- Consider an arbitrary sensor 
Sh € S' . Since St is the sensor in 5i2 with the smallest right extension in Cq, the right 
extension of Sh is larger than that of st in Cq. Depending on whether the covering intervals 
of St and Sh overlap in C, there are two subcases. 

If the covering intervals of st and Sh overlap in C, then let [0, x"] be the left-aligned interval 
that is covered by st and Sh in C (see Fig. E]). If we switch their relative positions by moving 
St to yt and moving y^ to x" — (i.e., the left extension of st is at and the right extension 
of Sh is at x"), then these two sensors still cover [0, x"] (see Fig. [3|), and thus the sensors in 5^' 
still form a coverage of [0, x"]. Further, after the above switch operation, the displacements of 
these two sensors are no bigger than A. To see this, clearly, the displacement of st is at most 
A. For the sensor Sh, it is easy to see that the switch operation moves Sh to the right. Since 
St covers p{x") in C, x" is no larger than the right extension of st in Cq, which is smaller than 
that of Sh in Cq. Thus, x" is smaller than x'^ + r/^, implying that the position of Sh after the 
switch operation is still to the left of its position in Cq. Hence, after the switch operation, the 
displacement of Sh is no bigger than A. In summary, after the switch operation, we obtain a 
new configuration C( in which the interval [0, x'] is still covered by the sensors of S[ and the 
position of the sensor st in is yt- 

If the covering intervals of St and Sh do not overlap in C, then suppose the sensors in the 
critical set Sc between Sh and st are s^, ■sj(i), -8/(2)1 • • • 1 •s/(m)i -^ii in the cover order. Clearly, 
the covering intervals of any two consecutive sensors in this list overlap in C. Below we show 
that we can switch the relative positions of st and such that we still form a coverage 

of [0,x'], and then we continue this switch procedure until st is switched with Sh- Note that 
since 5ii = 0, the right extension of for any 1 < j < m is larger than that of st in Cq. 

Let x'l be the maximum of and the left extension of in C , and X2 be the minimum of 
x' and the right extension of st in C (see Fig. [3]). Clearly, x'{ < X2. Thus, the sub-interval of 
[0, x'] covered by st and s/(m) in C is [x'{, X2]. We perform a switch operation on st and sj(m) 
by moving st to the left and moving S/(m) to the right such that the left extension of st is at 
x'l and the right extension of s/(m) is at X2 (see Fig. [3|). It is easy to see that after this switch 
operation, the sensors in Sc still form a coverage of [0,x']. Since the right extension of 
is larger than that of st in Cq, by a similar argument as above, we can also prove that after 
this switch, the displacements of both st and are no bigger than A. Then, we continue 

this switch process on st and sj(m_2), • • •, until st is switched with Sh, after which st 
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Figure 4: Illustrating the switch between st and s/{m)- 

is at ut, and we obtain a new configuration C[ in which the interval [0,x'] is still covered by 
the sensors in Sc ^ 5*^' and the position of the sensor st in C[ is y^. 

This completes the proof of the base case, i.e., there is always a configuration C[ in which the 
interval [0,x'] is covered by the sensors of S'^ and the position of the sensor s^^i) in €[ is yg(i)- 

We assume inductively that the claim holds for each k — 1 with 2 < k < i, i.e., there is a 
configuration C^_]^ in which the interval [0, x'] is covered by the sensors of S"- and the position of 
the sensor Sgi^j-j for each 1 < j < A; — 1 in C^„]^ is Ug^j)- In the following, we show that the claim 
holds for k, i.e., there is a configuration in which the interval [0, x'] is still covered by the sensors 
of S'^ and the position of the sensor Sg(j) for each 1 < j < /c in is yg[j)- The proof is quite similar 
to that for the base case and we only discuss it briefly below. 

Let t = g{k). If the position of st in C^_x is yt, then we are done (with = C^_x)- Otherwise, 
let y[ be the position of st in C'j^_^, with y[ ^ yt. Depending on st € Ski or st € Sk2i there are two 
cases. 

• If St € Ski, then yt = x[. Since yt is the rightmost position to which st is allowed to move 
and y'ti^yt, we have y't < yt- Depending on whether st is in the critical set Sc, there further 
are two subcases. 

If St ^ Sc, then the sensors in Sc always form a coverage of [0, x'] regardless of where st is. 
Thus, if we move st to yt, we obtain a new configuration from C^_i in which the sensors 
of S^ still form a coverage of [0, x'] and the position of the sensor s^q) for each 1 < j < fc in 
C'k is yg{j)- 

If St € Sc, then since yt > y't, if we move st from y[ to y^, is moved to the right. By Lemma 
djc), the interval [0,Rk-i] is covered by the sensors of ^^-i = {^^(i), Sg(2)) • • • > ^^(/c-i)} in 
C'^_^ (since they are in positions yg{i),yg{2), ■ ■ ■ ,yg{k-i), respectively). When st is at yt, st 
still covers the point p~^{Rk-i). Thus, after moving st to yt, we obtain a new configuration 
from C^_i in which the sensors of 5^' still form a coverage of [0, x']. 

• If St € Sk2, then S^i = in this case, and st is the sensor in Sk2 with the smallest right 
extension. If st ^ Sc, then by the same argument as above, we can obtain a configuration 
from C'f^_^ in which the interval [0, x'] is still covered by the sensors of S"- and the position of 
the sensor Sg(j) for each 1 < j < /c in is yg(j)- Below, we discuss the case when st S Sc- 

In 5c some sensors must cover the point p~^{Rk-i) in C. Let S' be the set of sensors in Sc 
that cover p+(i2fc_i) in C. If st £ S', then y[ < yt since yt is the rightmost position for st 
to cover In this case, again, by the same argument as above, we can move st to 

the right from y't to yt to obtain a configuration from C^_i in which the interval [0,x'] is 
still covered by the sensors of S'^. Otherwise (i.e., st ^ S'), consider a sensor Sh in S'. Let 
the sensors in Sc between Sh and st in the cover order be s^, sj(i), Sj-(2), . . . , sj(m), sj (this 
sequence may contain only Sh and st). Note that for each 1 < j < fc — 1, the sensor Sg(^j) is 
not in this sequence. Then by using a similar sequence of switch operations as for the base 
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case, we can obtain a new configuration from C^,_i such that the sensors of S'^ still form 
a coverage of [0, x']. Again, the position of the sensor for each 1 < j < /c in is yg(j)- 

This proves that the claim holds for k. Therefore, the claim is true. The lemma can then be 
easily proved by using this claim, as follows. 

Suppose the largest left-aligned interval that can be covered by the sensors of is [0, x']. Then 
by the above claim, there always exists a configuration C* for S'^ in which the interval [0, x'] is also 
covered by the sensors of S'^ and for each 1 < j < i, the position of the sensor in C* is yg(j)- 
Recall that Si = {sg(i), ■5^(2) i • • • i Then for each sensor st € S'- \ Si, the rightmost point that 

can be covered by st is x'^ + rt- Recall that in the configuration Cj, for each 1 < j < i, the position 
of the sensor s^q) is yg(j), and for each sensor st € S^ \ Si, the position of st is x[. Further, by the 
definition of S'^, the right extensions of all sensors in S'^ are at most Ri in Cj. Therefore, the right 
extensions of all sensors in S'^ are also at most Ri in C*, implying that x' < Ri. On the other hand, 
by Lemma mc), the sensors of Si form a coverage of [0, in C* . Thus, [0,x'] = [0,Ri], and the 
lemma follows. □ 

Finally, we prove the correctness of our algorithm based on Lemma [2j Suppose our algorithm 
reports A < A* in step i. Then according to the algorithm, Ri-i < L and both Sn and Si2 are 0. 
Let Sl_^ be the set of sensors whose right extensions are at most Ri-i in Cj_i. Since both Sn and 
Si2 are 0, no sensor in S\Sl_-^ can cover any point to the left of the point (and including 

By Lemma [21 [0, is the largest left-aligned interval that can be covered by the 

sensors of 5'j'_^. Hence, the sensors in S cannot cover the interval [0,p^(i?i_i)]. Due to Ri-i < L, 
we have [0,p^(i?j_i)] C [0,L]; thus the sensors of S cannot cover B = [0,L]. In other words, there 
is no feasible solution for the distance A. This establishes the correctness of our algorithm. 

2.1.4 The Algorithm Implementation 

For the implementation of the algorithm, we first discuss a straightforward approach that runs in 
0(n log n) time. Later, we give another approach which, after 0(n log n) time preprocessing, can 
determine whether A* < A in 0(n) time for any given A. Although the second approach does not 
change the overall running time of our decision algorithm, it does help our optimization algorithm 
in Section 12.21 to run faster. 

In the beginning, we sort the 2n extensions of all sensors by the x-coordinate, and move each 
sensor Sj G 5 to x'- to produce the initial configuration Cq. During the algorithm, for each step i, we 
maintain two sets of sensors, Sn and Si2, as defined earlier. To this end, we sweep along the x-axis 
and maintain Sn and Si2, using two sweeping points pi and p2, respectively. Specifically, the point pi 
follows the positions Rq {= 0), Ri, R2, ■ ■ ., and p2 follows the positions Ro + 2X, Ri + 2X, R2 + 2X, .... 
Thus, p2 is kept always by a distance of 2A to the right of pi. To maintain the set Sn, when the 
sweeping point pi encounters the left extension of a sensor, we insert the sensor into Sn; when pi 
encounters the right extension of a sensor, we delete the sensor from Sn. In this way, when the 
sweeping point pi is at Ri-i, we have the set Sn ready. To maintain 5,2, the situation is slightly 
more subtle. First, whenever the sweeping point p2 encounters the left extension of a sensor, we 
insert the sensor into Si2- The subtle part is at the deletion operation. By the definition of Si2, if 
the left extension of any sensor is less than or equal to Ri-i, then it should not be in Si2. Since 
eventually the first sweeping point pi is at Ri-\ in step i, whenever a sensor is inserted into the 
first set Sn, we need to delete that sensor from Si2 (note that this sensor must be in Si2 since p2 
is to the right of pi). Thus, a deletion on Si2 happens only when the same sensor is inserted into 
Sii. In addition, we need a search operation on Sn for finding the sensor in Sn with the largest 
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Figure 5: Illustrating the search for the sensor s: the sensors before pi are removed from Sr; the two sensors below 
s are redundant sensors for this step. 

right extension, and a search operation on Si2 for finding the sensor in Si2 with the smahest right 
extension. 

It is easy to see that there are 0{n) insertions and deletions in the entire algorithm. Further, 
the search operations on both Sn and Si2 are dependent on the right extensions of the senors. By 
using a balanced binary search tree to represent each of these two sets in which the right extensions 
of the sensors are used as keys, the algorithm runs in 0(n log n) time. 

In the sequel, we give the second approach which, after O(nlogn) time preprocessing, can 
determine whether A* < A in 0{n) time for any given A. 

In the preprocessing, we compute two sorted lists Sl and Sr, where Sl contains all sensors 
sorted by the increasing values of their left extensions and Sr contains all sensors sorted by the 
increasing values of their right extensions. Consider any value A. Later in the algorithm, for each 
step i, our algorithm will determine the sensor s^^j) by scanning the two lists. We will show that 
when the algorithm finishes, each sensor in Sl is scanned at most once and each sensor in Sr is 
scanned at most three times, and therefore, the algorithm runs in 0(n) time. 

Initially, we move each sensor Si £ S to x[ to produce the initial configuration Cq. During the 
algorithm, we sweep along the x-axis, using a sweeping point pi. Specifically, the point pi follows 
the positions Rq (=0), Ri, R2, ■ ■ ■■ With a little abuse of notation, we also let pi be the coordinate 
of the current position of pi. Initially, pi = 0. 

Consider a general step i and we need to determine the sensor Sg(i)- In the beginning of this 
step, pi is at the position Ri-i. We scan the list Sl from the beginning until the left extension 
of the next sensor is strictly to the right of pi . For each scanned sensor Sj , if its right extension 
is strictly to the right of pi, then it is in Sn by the definition of Sn. Thus, the above scanning 
procedure can determine Sn, after which we can easily find the sensor 5^(1) in Sn if Sn ^ 0. In fact, 
we can compute Sg(i) directly in the above scanning procedure. In addition, for any sensor in Sl 
that is scanned, we remove it from Sl (and thus it will never be scanned later in the algorithm any 
more). If Sn 7^ 0, then Sg(j) is determined and we move the sweeping point pi to the right extension 
of (i.e., pi = Ri). If pi > L, we terminate the algorithm and report A* < A; otherwise, we 
continue on the step i + 1. Below, we discuss the case Sn = 0. 

If Sii = 0, then s^^j) is the sensor in 5^2 with the smallest right extension if 5,2 7^ 0- Specifically, 
among the sensors (if any) whose left extensions are larger than pi (= Ri-i) and at most pi + 2A, 
Sg{i) is the sensor with the smallest right extension. To find s^^j), we scan the list Sr from the 
beginning until we find the first sensor s whose left extension is larger than pi and at most pi + 2A 
(see Fig. [5]). If such a sensor s does not exist, then 5^2 = 0, and we terminate the algorithm and 
report A* > A. Below, we assume we have found such a sensor s. Since the sensors in Sr are sorted 
by their right extensions, Sg(j) is exactly the sensor s. Further, unlike the scanning on Sl where 
each scanned sensor is removed immediately, for each scanned sensor in Sr, we remove it only if 
its right extension is to the left of pi (see Fig. [5]). It is easy to see that the removed sensors (if 
any) are consecutive from the beginning of Sr. Let Sr be the list after all removals. If (= s) 
is not the first sensor in Sr, then for any sensor Sj in Sr before the left extension of Sj must 
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be larger than pi + 2A; we call the sensors in Sr before Sg(j) the redundant sensors for the step 
i (see Fig. [5]). Later we will show that these sensors will be not redundant any more in the later 
algorithm. In summary, for each sensor scanned in the original Sr in this step, it is removed, or 
a redundant sensor, or Finally, we move to the left such that its left extension is at pi, 

and then we move pi to the right extension of Sg(j) (i.e., pi = Ri). If pi > L, we terminate the 
algorithm and report A* < A; otherwise, we continue on the next step i + 1. 

To analyze the algorithm, it is easy to see each sensor in this list Sl is scanned at most once. 
For the list Sn, this may not be the case as the redundant sensors may be scanned again in the 
later algorithm. However, the following lemma shows that this would not be an issue. 

Lemma 3 If a sensor Sj is a redundant sensor for the step i, then it will be not a redundant sensor 
again in the later algorithm. 

Proof: Consider the moment right after the step i. The sweeping point pi is at the right extension 
of Sg(i)- To prove the lemma, since pi always moves to the right, by the definition of the redundant 
sensors, it is sufficient to show that the left extension of Sj is at most pi + 2A, as follows. 

Consider the moment in the beginning of the step i (the sensor has not been moved to the 
left). Since Sj is a redundant sensor for the step i, the sensor Sg(j) is from Si2 and the left extension 
of Sg(j) is at most Ri-i + 2A. Thus, the right extension of is at most Ri-i + 2rg(j) + 2A. Recall 
that the right extension of sj is less than that of (since Sj is before in Sr). Therefore, the 
right extension of sj is at most Ri-i + '^'''g{i) + 2A. Now consider the moment right after the step i. 
The sweeping point pi is at the position Ri-i + '^fg{i)- Hence, the right extension of Sj is at most 
pi + 2A, which implies that the left extension of Sj is at most pi + 2A. The lemma thus follows. □ 

The preceding lemma implies that any sensor can be a redundant sensor in at most one step. 
Therefore, for the list Sr, each sensor is scanned at most once when it is removed, once as a 
redundant sensor, and once when it is found as Thus, each sensor in Sr is scanned at most 

three times. Hence, after the two lists Sl and Sr are obtained, the running time of the algorithm 
is 0{n). 

Theorem 1 After 0(n log n) time preprocessing, for any A, we can determine whether X* < X in 
0{n) time; further, if X* < A, we can compute a feasible solution in 0{n) time. 

2.1.5 Another Decision Version 

Our optimization algorithm in Section [2.21 also needs to determine whether A* is strictly less than 
A (i.e.. A* < A) for any A. By modifying our algorithm for Theorem[Tl we have the following result. 

Theorem 2 After 0{n\ogn) time preprocessing, for any value X, we can determine whether X* < X 
in 0{n) time. 

Proof: We first apply the algorithm for Theorem[T]on the value A. If the algorithm reports A* > A, 
then we know A* < A is false. Otherwise, we have A* < A. In the following, we modify the algorithm 
for Theorem [T] to determine whether A* < A, i.e.. A* is strictly smaller than A. Note that this is 
equivalent to deciding whether A* < A — e for any arbitrarily small constant e > 0. Of course, we 
cannot enumerate all such small values e. Instead, we add a new mechanism to the algorithm for 
Theorem [1] such that the resulting displacement of each sensor is strictly smaller than A. 

At the start of the algorithm, we move all sensors to the right by a distance A to obtain the 
configuration Cq. But, the displacement of each sensor should be strictly less than A. To ensure 
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this, later in the algorithm, if the destination of a sensor Sj is set as yi = x'-, then we adjust this 
destination of Si by moving it to the left slightly such that Sj's displacement is strictly less than A. 

Consider a general step i of the algorithm. We define the set Su in the same way as before, i.e., 
it consists of all sensors covering the point p'^{Ri-i) in Q-i. If Sn 7^ 0, then the algorithm is the 
same as before. In this case, the sensor s^(j-) chosen in this step has a displacement of exactly A, 
which is actually "illegal" since the displacement of each sensor should be strictly less than A. We 
will address this issue later. However, if Sn = 0, then the set Si2 is defined slightly different from 
before. Here, since Sn = 0, we have to use a sensor to the right of Ri-i in Ci-i to cover 
Since the displacement of each sensor should be strictly less than A, we do not allow any sensor 
to move to the left by exactly the distance 2A. To reflect this difference, we define Si2 as the set 
of sensors in Cj-i each of which has its left extension larger than Ri-i and strictly smaller than 
Ri-i + 2 A (previously, it was "at most"). In this way, if we move a sensor in Si2 to the left to cover 
then the displacement of that sensor is strictly less than A. The rest of the algorithm is 
the same as before. We define the type I and type II sensors in the same way as before. 

If the algorithm terminates without finding a feasible solution, then it must be A* > A; otherwise, 
the algorithm finds a "feasible" solution SOL with a critical set S'^ = ■Sg(2)) • • • > Sg{m)}- I^^t, 

this does not necessarily mean A* < A since in SOL, the displacements of some sensors in S*^ may 
be exactly A. Specifically, all type I sensors in S^ are in the same positions as they are in Cq and 
thus their displacements are exactly A. In contrast, during the algorithm, the type II sensors in S^ 
have been moved strictly to the left with respect to their positions in Cq; further, due to our new 
definition of the set Si2, the displacements of all type II sensors are strictly less than A. Therefore, 
if there is no type I sensor in S'^, then the displacement of each sensor in S'^ is strictly less than A 
and thus we have A* < A. Below we assume S**^ contains at least one type I sensor. To make sure 
that A* < A holds, we need to find a real feasible solution in which the displacement of each sensor 
in S is strictly less than A. On the other hand, to make sure that A* > A holds, we must show that 
there is no real feasible solution. For this, we apply the following algorithmic procedure. 

We seek to adjust the solution SOL to produce a real feasible solution. According to our 
algorithm, for each sensor Sj € 6"^, if it is a type I sensor, then yi = x'- and thus its displacement 
is exactly A; otherwise, its displacement is less than A. The purpose of our adjustment of SOL is 
to move all type I sensors slightly to the left so that (1) their displacements are strictly less than 
A, and (2) we can still form a coverage of B. In certain cases, we may need to use some sensors in 
S* \ S"*^ as well. Also, we may end up with the conclusion that no real feasible solution exists. 

According to our algorithm, after finding the last sensor ^^(to) in S'^, we have Rm > L. If 
Rm > L, then we can always adjust SOL to obtain a real feasible solution by shifting each sensor 
in to the left by a very small value e such that (1) the resulting displacement of each sensor in 
iS^ is less than A, and (2) the sensors of S"^ still form a coverage of B. Note that there always exists 
such a small value e such that the above adjustment is possible. Therefore, if Rm > L, then we 
have A* < A. 

If Rm = L, however, then the above strategy does not work. There are two cases. If there is a 
sensor st € S \ S'^ such that xt G {L — X — rt, L + X + rt), then wc can also obtain a real feasible 
solution by shifting the sensors of S'^ slightly to the left as above and using the sensor st to cover 
the remaining part of B around L that is no longer covered by the shifted sensors of S'^; thus we 
also have A* < A. Otherwise, we claim that it must be A* > A. Below we prove this claim. 

Consider the rightmost Type I sensor Sj in S"^. Suppose Si = i.e., Si is determined in step 
j. Thus, Si is at x[ in SOL. Let e > be an arbitrarily small value (we will determine below 
how small it should be). Since we have assumed that the extensions of all sensors are different, the 
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value e can be made small enough such that by moving Sj to — e in Co, the relative order of the 
extensions of all sensors remains the same as before. Further, according to our algorithm above, 
the value e can also be small enough such that the behavior of the algorithm is the same as before, 
i.e., the algorithm finds the same critical set with the same cover order as before. It is easy to 
see that such a small value e always exists. Note that our task here is to prove our claim A* > A is 
true, and knowing that such a value e exists is sufficient for our purpose and we need not actually 
find such a value e in our algorithm. 

Now, in step j, the new value Rj, which is the right extension of Sj, is e smaller than its value 
before since Si was at x'^ in Cq. Because Si is the rightmost type I sensor in S'^, after step j, all 
sensors in 5"^ determined after Si (if any) are of type II and thus are moved to the left such that they 
are all in attached positions along with s,, which implies that the right extension of the last sensor 
Sg(^rn) ill S'^ is s-lso € Smaller than its previous value (which was L). Hence, after step m, the sensors 
in Sm covers [0, L — e]. As discussed above, if e is made small enough, the behavior of the algorithm 
is the same as before. By a similar analysis, we can also establish a result similar to Lemma [2j 
Namely, [0, L — e] is the largest left-aligned interval that can be covered by the sensors in S'^ in this 
setting (here, S'^^ is the set of sensors whose right extensions are at most L — e in the configuration 
after step m). We omit the detailed analysis for this, which is very similar to that for Lemma [2j 
Note that S*^ = Sm- Since there is no sensor st & S\S'^ such that xt € {L — X — rt, L + X + rt), the 
interval (L — e, L] cannot be fully covered by the sensors in S. The above discussion implies that 
if we do not allow the displacement of Si to be larger than A — e, then there would be no feasible 
solution even if we allow the displacements of some other sensors (i.e., those type I sensors in S'^ 
before Si, if any) to be larger than A — e (but at most A). Thus, A* < A — e cannot be true. That is, 
A* > A — e holds. Further, it is easy to see that, by a similar argument, for any fixed value e' > 
with e' < e, we also have A* > A — e'. Hence, we obtain A* > A. 

This finishes the discussion on how to determine whether A* < A. It is easy to see that the 
above algorithm can also be implemented in 0(n) time for each value A, after 0(n log n) time 
preprocessing. The theorem thus follows. □ 

Theorems [1] and [2] together lead to the following corollary. 

Corollary 1 After O(nlogn) time preprocessing, for any value X, we can determine whether X* = 
A in 0{n) time. 

2.2 The Optimization Version of the General BCLS 

In this section, we discuss the optimization version of the general BCLS problem. We show that it 
is solvable in 0(n^ log n) time, thus settling the open problem in [8]. 

The main difficulty is that we do not know the order of the critical sensors in an optimal 
solution. Our strategy is to determine a critical set of sensors and their cover order in a feasible 
solution for the (unknown) optimal value A*. The idea is somewhat similar to parametric search 
[6|ll4j and here we "parameterize" our algorithm for Theorem [TJ But, unlike the typical parametric 
search [6l |T3j, our approach does not involve any parallel scheme and thus is practical. We first 
give an overview of this algorithm. In the following discussion, the "decision algorithm" refers to 
our algorithm for Theorem [1] unless otherwise stated. 

Recall that given any value A, step i of our decision algorithm determines the sensor and 
obtains the set Si = {•Sg(i), ■5g(2), . . . , ■Sg(j)}, in this order, which we also call the cover order of the 
sensors in Si. In our optimization algorithm, we often use A as a variable. Thus, Si{X) (resp., Ri(X), 
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Sg(j)(A), and Ci(A)) refers to the corresponding Si (resp., Ri, and Q) obtained by running 

our decision algorithm on the specific value A. Denote by Cj the configuration of the input. 

Our optimization algorithm takes at most n steps. Step i receives an interval Xf_^) and 

a sensor set 5j_i(A*), with the algorithm invariants that A* € {^j-i, ^i-i) (although we do not 
know the value A*) and for any value A G (Aj^^^, A?_^), we have 5j_i(A) = 5j_i(A*) and their 
cover orders are the same. Step i either finds the value A* or determines a sensor Sg(j)(A*). The 
interval (Aj_^, Xf_^) will shrink to a new interval {Xj,Xf) C iXj_^,\f_^) and we also obtain the set 
Si{X*) = Si-i{X*) U {■Sg(j)(A*)}. Each step can be performed in 0(n log n) time. The details of the 
algorithm are given below. 

Initially, let 5o(A*) = 0, i?o(A*) = 0, Aj = 0, and A^ = +oo. 

Consider a general step i for i > 1 and we have the interval {Xl_i, Xj_i) and the set 5j_i(A*). 
While discussing the algorithm, we will also prove inductively the following lemma about the 
function Ri{X) with variable A G {Xj,Xf). 

Lemma 4 (a) The function Ri{X) for X G (A]^,A^) is a line segment of slope 1 or 0. (b) We 
can compute the function Ri{X) for X G {Xj,Xf) explicitly in 0{n) time, (c) Ri{X) < L for any 
AG(Ai,A2). 

In the base case for i = 0, the statement of Lemma H] obviously holds. We assume the lemma 
statement holds for i — 1. We will show that after step i with i > 1, the lemma statement holds for 
i, and thus the lemma will be proved. 

Again, in step i, we need to determine the sensor Sg(j)(A*) and let 5j(A*) = 5j_i(A*)U{sg(j)(A*)}. 
We will also obtain an interval {Xj,Xf) such that A* G (Aj^,A?) C {Xj_i,Xi_i) and for any A G 
(AJ,A?), Si{X) = Si{X*) holds (with the same cover order). The details are given below. We 
assume that we already compute explicitly the function i?i_i(A) for A G {Xj_^, X^_^), which takes 
0{n) time by our assumption that the statement of Lemma S] holds for i — 1. 

To find the sensor Sg(j)(A*), we first determine the set 5ji(A*). Recall that Sii{X*) consists 
of all sensors covering the point p~^ {Ri-i{X*)) in the configuration Cj_i(A*). For each sensor in 
S \ Si-i{X*), its position in the configuration Ci_i(A) with respect to A G [Xj^^, Xf_^) is a function 
of slope 1. As A increases in {Xj_i, Xf_i), by our assumption that LemmalD^a) holds for i — 1, the 
function i?j_i(A) is a line segment of slope 1 or 0. If i2i_i(A) is of slope 1, then the relative position 
of i?j_i(A) in Cj_i(A) does not change and thus the set Sii{X) does not change; if the function 
Ri-i{X) is of slope 0, then the relative position of i?j„i(A) in Cj_i(A) is monotonically moving to 
the left. Hence, there are 0{n) values for A in (A^^^^, Xf_-^^) that can incur some changes to the set 
Sii{X) and each such value corresponds to a sensor extension; further, these values can be easily 
determined in O(nlogn) time by a simple sweeping process (we omit the discussion of it). Let 
Aji be the set of all these A values. Let Aji also contain both Xj_^ and Xf_^, and thus, Xj_^ and 
^f-i are the smallest and largest values in Aji, respectively. We sort the values in An. For any 
two consecutive values Ai < A2 in the sorted Aji, the set Sii{X) for any A G (Ai,A2) is the same. 
By using binary search on the sorted Aji and our decision algorithm, we determine (in O(nlogn) 
time) the two consecutive values Ai and A2 in Aji such that Ai < A* < A2. Further, by Corollary [H 
we determine whether A* = A2. If A* = A2, then we are done. Otherwise, based on our discussion 
above, Sji(A*) = 5'ji(A) for any A G (Ai,A2). Thus, to compute 5ji(A*), we can pick an arbitrary 
A in (Ai,A2) and find 5'ii(A) in the same way as in our decision algorithm. Hence, 5ji(A*) can be 
easily found in O(nlogn) time. Note that A* G (Ai,A2) C {Xj_^, Xf_^). 

If Sii{X*) 7^ 0, then s^(j-)(A*) is the sensor in 5ii(A*) with the largest right extension. An obvious 
observation is that for any A G (Ai, A2), the sensor in Sii{X*) with the largest right extension is the 
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same, which can be easily found. We let Xj = Ai and Xj = X2- Let Si{X*) = Si-i{X*) U {sg(i)(A*)}. 
The algorithm invariants hold. Further, as A increases in (Ai,A2), the right extension of Sg(j)(A), 
which is Ri{X), increases by the same amount. That is, the function Ri{X) on (Ai,A2) is a line 
segment of slope 1. Therefore, we can compute Ri{X) on (Ai,A2) explicitly in constant time. This 
also shows Lemma m^a) and (b) hold for i. 

Because the function Ri{X) on (Ai,A2) is a line segment of slope 1, there are three cases de- 
pending on the values Ri{X) and L: (1) Ri{X) < L for any A G (Ai,A2), (2) Ri{X) > L for any 
A € (Ai,A2), and (3) there exists a unique value A' € (Ai,A2) such that Ri{X') = L. For Case (1), 
we proceed to the next step, along with the interval {Xj, A?). Clearly, the algorithm invariants hold 
and Lemma m^c) holds for i. For Case (2), the next lemma shows that it actually cannot happen 
due to A* G (Ai, A2). 

Lemma 5 It is not possible that Ri{X) > L for any X G (Ai, A2). 

Proof: Assume to the contrary that Ri{X) > L for any A G (Ai, A2). Since A* G (Ai, A2), let A" be 
any value in (Ai,A*). Due to A" G (Ai,A2), we have Ri{X") > L. But this would implies that we 
have found a feasible solution where the displacement of each sensor is at most A", which is smaller 
than A*, incurring contradiction. □ 

For the Case (3), since Ri{X) on (Ai,A2) is a line segment of slope 1, we can determine in 
constant time the unique value A' G (Ai, A2) such that Ri{X') = L. Clearly, A* < A'. By Corollary 
[H we determine whether A* = A'. If A* = A', then we are done; otherwise, we have A* G (Ai,A') 
and update A? to A'. We proceed to the next step, along with the interval {Xj,Xf). Again, the 
algorithm invariants hold and Lemma [ll^c) holds for i. 

If Sii{X*) = 0, then we need to compute Si2{X*). For any A G (Ai,A2), the set Si2{X) consists 
of all sensors whose left extensions are larger than i?.j_i(A) and at most i?i_i(A) + 2A in the 
configuration Cj_i(A). Recall that the function i?j_i(A) on {Xl_i, Xf_i) is linear with slope 1 or 
0. Due to (Ai,A2) C {Xl_i, Xf^^), the linear function i?i_i(A) + 2A on (Ai,A2) is of slope 3 or 2. 
Again, as A increases, the position of each sensor in 5 \ 5j_i(A*) in Ci_i(A) is a linear function of 
slope 1. Therefore, there are 0(n) A values in (Ai, A2) each of which incurs some change to the set 
5*42 (A) and each such A value corresponds to a sensor extension. Further, these values can be easily 
determined in 0(n log n) time by a sweeping process (we omit the discussion for this). (Actually, 
as A increases, the size of the set 5'j2(A) is monotonically increasing.) Let Aj2 denote the set of 
these A values, and let Aj2 contain both Ai and A2. Again, |Aj2| = 0{n). We sort the values in 
A,2. Using binary search on the sorted Aj2 and our decision algorithm, we determine (in 0(n log n) 
time) the two consecutive values A'^^ and A2 in Aj2 such that X'l < X* < A2. Further, by Corollary 
[H we determine whether A* = A2. If A* = A2, then we are done. Otherwise, Si2{X*) = Si2{X) for 
any A G (A'j^jAg), which can be easily found. Note that A* G (A'^,A2) C (Ai,A2). 

After obtaining 5^2 (A*), Sg(j)(A*) is the sensor in S'i2(A*) with the smallest right extension. As 
before, the sensor in Si2{X) with the smallest right extension is the same for any A G (A'^, A2). Thus, 
Sg(j)(A*) can be easily determined. We let Xj = X'^ and Xf = A2. Let Si{X*) = S'j_i(A*)U{s£,(j)(A*)}. 
The algorithm invariants hold. Further, we examine the function Ri{X), i.e., the right extension 
of Sc,(j)(A) in the configuration Cj(A), as A increases in (A'j^, A2). Since Sg(j„i)(A*) and Sg(j)(A*) are 
always in attached positions in this case, for any A G (A']^,A2), we have Ri{X) = Ri-i{X) + 2rg(^iy 
Thus, the function Ri{X) is a vertical shift of i?i_i(A) by the distance "^rg^^^y Because we already 
know explicitly the function i?j_i(A) for A G {X'i,X'2), which is a line segment of slope 1 or 0, the 
function Ri{X) can be computed in constant time, which is also a line segment of slope 1 or 0. Note 
that this shows that Lemma |3)^a) and (b) hold for i. 
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Similarly to the case when Sii{X*) ^ 0, since the function Ri{X) in (A'^,A2) is a line segment 
of slope 1 or 0, there are three cases depending on the values Ri{X) and L: (1) Ri{X) < L for any 
A G {X'i,X2), (2) Ri{X) > L for any A G {Xi,X2), and (3) there exists a unique value A" G (A'i,A2) 
such that Ri{X") = L. For Case (1), we proceed to the next step, along with the interval {X\^X^). 
Clearly, the algorithm invariants hold and Lemma Ul^c) holds for i. Similarly to Lemma [5l Case (2) 
cannot happen due to A* G (A'^,A2). For the Case (3), since Ri{X) in {X'i,X'2) is a line segment of 
slope 1 or 0, we can compute in constant time the unique value A" G (A'^, Ag) such that Ri{X") = L. 
Clearly, A* < A". By Corollary [U we determine whether A* = A". If A* = A", we are done; 
otherwise, we have A* G (A'^, A") and update A? to A". We proceed to the next step, along with the 
interval (AJ,A|). Again, the algorithm invariants and Lemma|3|^c) hold for i. 

This finishes the discussion of step i of our algorithm. Note that in each case where we proceed 
to the next step. Lemma U] holds for i, and thus Lemma [5] has been proved. The running time of 
step i is clearly bounded by 0(n log n). 

In at most n steps, the algorithm will stop and find the value A*. Then by applying our decision 
algorithm on A = A*, we finally produce an optimal solution in which the displacement of every 
sensor is at most A*. Since each step takes O(nlogn) time, the total time of the algorithm is 
0(n^ logn). 

Theorem 3 The general BCLS problem is solvable in O(n^logn) time. 

We shall make a technical remark. The typical parametric search [SI E] usually returns with 
an interval containing the optimal value and then uses an additional step to find the optimal value. 
In contrast, our algorithm is guaranteed to find the optimal value A* directly. This is due to the 
mechanism in our algorithm that requires Ri{X) < L for any A G (A[,A?) after each step i if the 
algorithm is not terminated. This mechanism actually plays the role of the additional step used in 
the typical parametric search. 

3 The Uniform Case of BCLS 

In this section, we present an O(nlogn) time algorithm for the uniform case of BCLS. Previously, 
the best known algorithm for it takes O(n^) time [8]. Further, for the special uniform case when 
all sensors are initially located on the barrier B, we solve it in 0(n) time. 

3.1 Preliminaries 

Recall that in the input, all sensors are ordered from left to right by their initial positions, i.e., 
xi < X2 < ■ ■ ■ < Xn- Suppose in a solution, the destination of each sensor Sj is (1 < i < ra); then 
we say that the solution is order preserving if yi < y2 ^ • • • ^ Vn- In the uniform case, since all 
sensors have the same range, we let r denote the sensor range. The next lemma was known [8]. 

Lemma 6 (Czyzowicz et al. [8\) For the uniform case, there is always an optimal solution that is 
order preserving. 

As discussed in [8], Lemma[6]is not applicable to the general BCLS. Consequently, the approach 
in this section does not work for the general BCLS. 

Based on the order preserving property in Lemma [6l the previous O(n^) time algorithm [8] tries 
to cover B from left to right; each step picks the next sensor and re-balances the current maximum 
sensor movement. Here, we take a completely different approach. 
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Denote by A* the maximum sensor movement in an optimal solution. We use OPT to denote 
an optimal order preserving solution in which the destination for each sensor Si is yi {1 < i < n). 
For each sensor Sj, if Xi > yi (resp., Xi < yi), then we say Sj is moved to the left (resp., right) by a 
distance \xi — yi\. A set of sensors is said to be in attached positions if the union of their covering 
intervals is a continuous interval on the x-axis whose length is equal to the sum of the lengths of 
the covering intervals of these sensors. A single sensor is always in attached position. The following 
lemma was proved in 

Lemma 7 (Czyzowicz et al. fH]) 7/ A* > 0, then in OPT, there exist a sequence of consecutive 
sensors Si, Sj+i, . . . ,Sj with i < j such that they are in attached positions and one of the following 
three cases is true, (a) The sensor sj is moved to the left by the distance A* and yi = r (i.e., the 
sensors Sj, Sj+i, . . . , Sj together cover exactly the interval [0, 2r(j — i + 1)]). (b) The sensor Si is 
moved to the right by the distance A* and yj = L — rj. (c) For i ^ j (i.e., i < j), the sensor Si is 
moved to the right by the distance A* and the sensor Sj is moved to the left by the distance A*. 

Cases (a) and (b) in Lemma [7] are symmetric. By Lemma [71 for each pair of sensors Sj and sj 
with i < j, we can compute three distances Ai(i, j), A2(i, j), and A3(i, j) corresponding to the three 
cases in Lemma[7]as candidates for the optimal distance A*. Specifically, Xi{i,j) = Xj — [2r{j—i)+r], 
where the value 2r(j — z) + r is supposed to be the destination of the sensor Sj in OPT if case (a) 
holds. Symmetrically, \2{i-,j) = [L — 2r(j — i) — r] — Xj. Let Xsii^j) = [xj — Xi — 2r(j — i)]/2 for 
i < j. Let A be the set of all Ai(i, j), A2(i,i), and X3{i,j) values. Clearly, A* S A and |A| = 0(n^). 
By using an algorithm for the decision version of the uniform case to search in A, one can find 
the value A*. Recall that the decision problem is that given any value A, determine whether there 
exists a feasible solution for covering B such that the moving distances of all sensors are at most 
A. Thus, A* is the smallest value in A such that the answer to the decision problem on that value 
is "yes". A simple greedy 0{n) time algorithm was given in |8] for the decision problem. 

Lemma 8 (Czyzowicz et al. [8j) The decision version of the uniform case is solvable in 0{n) time. 

But, the above approach would take Q.(n?) time due to |A| = Q(n?). To reduce the running 
time, we cannot compute the set A explicitly. In general, our 0(n log n) time algorithm uses the 
following idea. First, instead of computing all elements of A explicitly, we compute one element 
of A whenever we need it (we may do some preprocessing for this). Second, suppose we already 
know (implicitly) a sorted order of all values in A; then we can use binary search and the decision 
algorithm for Lemma [S] to find A*. However, we are not able to order the values of A into a single 
sorted list; instead, we order them (implicitly) in 0(n) sorted lists and each list has 0{n) values. 
Consequently, by a technique called binary search in sorted arrays [5], we compute A* in 0(n log n) 
time. The details of our algorithm are given in the next subsection. 

3.2 Our Algorithm for the Uniform Case 

Due to the order preserving property, it is easy to check whether A* = in 0(n) time. In the 
following, we assume A* > 0. 

We focus on how to order (implicitly) the elements of A into 0{n) sorted lists and each list 
contains 0{n) elements. We also show that after preprocessing, each element in any sorted list can 
be computed in 0(1) time using the index of the element. We aim to prove the next lemma. 
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Lemma 9 In 0(n log n) time, the elements of A can be ordered (implicitly) into 0{n) sorted lists 
such that each list contains 0{n) elements and each element in any list can be computed in constant 
time by giving the index of the list and the index of the element in the list. 

The following technique, called binary search on sorted arrays [5], will be applied. Suppose there 
is a "black-box" decision procedure 11 available such that given any value a, H can report whether 
a is a feasible value to a certain problem in 0{T) time, and further, if a is a feasible value, then any 
value larger than a is also feasible. Given a set of m arrays Ai, 1 < i < m, each containing 0{n) 
elements in sorted order, the goal is to find the smallest feasible value 6 in A = U^^^dj. Suppose 
each element of any array can be obtained in constant time by giving its indices. An algorithm for 
the following result was presented by Chen et al. in [5\. 

Lemma 10 (Chen et al. [5J) The value 6 in A can be computed in 0((m + T) log (nm)) time. 

If we use the algorithm for Lemma [H] as the decision procedure 11, then by Lemmas M and \TU\ 
we can find A* in A in 0(n log n) time. After A* is found, we can apply the algorithm for Lemma 
[8] to compute the destinations of all sensors, in 0{n) time. Hence, we have the following result. 

Theorem 4 The uniform case of the BCLS problem is solvable in O(nlogn) time. 

In the rest of this subsection, we focus on proving Lemma EJ 

For each 1 < t < 3, let At denote the set of all Xt{i,j) values. Clearly, A = Ai U A2 U A3. We 
seek to order each of the three sets Ai, A2, and A3 into sorted lists. 

We discuss Ai first. This case is trivial. It is easy to see that for each given value j, we 
have Ai(ii,j) < Ai(i2,j) for any ii < ^2 < J. Thus, for every value j, we have a sorted list 
Ai(l, j), Ai(2, j), . . . , Ai(j, j) of j elements, and each element can be computed in constant time by 
using the index of the element in the list. Therefore, we have n sorted lists, and clearly, the set of 
elements in all these lists is exactly Ai. Hence we have the following lemma. 

Lemma 11 In O(nlogn) time, the elements of Ai can be ordered (implicitly) into 0{n) sorted 
lists such that each list contains 0{n) elements and each element in any list can be computed in 
constant time by giving the index of the list and the index of the element in the list. 

The set A2 can be processed in a symmetric manner as Ai, and we omit the details. 

Lemma 12 In 0(n log n) time, the elements 0/ A2 can be ordered (implicitly) into 0{n) sorted 
lists such that each list contains 0{n) elements and each element in any list can be computed in 
constant time by giving the index of the list and the index of the element in the list. 

In the following, we focus on ordering (implicitly) the set A3 and showing the following lemma, 
which, together with Lemmas 1111 and 1121 proves Lemma [9l 

Lemma 13 In 0(n log n) time, the elements of A3 can be ordered (implicitly) into 0{n) sorted 
lists such that each list contains 0{n) elements and each element in any list can be computed in 
constant time by giving the index of the list and the index of the element in the list. 
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Figure 6: Illustrating three groups of sensors: Gi-i,Gi, and d+i. The sensors with indices ai-i,ai, and a^+i 
(resp., bi-i,bi, and fei+i) are in Sa (resp., Sb), which are indicated by thick segments. 

Proving Lemma [T3l is a main challenge to our uniform case algorithm. The reason is that, unlike 
Ai and A2, for a given j, for any I < h < 12 < j, either XsikJ) < hihj) or XsihJ) > XsihJ) 
is possible. Hence, to prove Lemma [131 we have to find another way to order the elements of A3. 

Our approach is to first remove some elements from A3 that are surely not A* (for example, 
negative values cannot be A*). We begin with some intuitions. We say two intervals on the x-axis 
are strictly overlapped if they contain more than one common point. In the following discussion, 
the sensors are always at their input positions unless otherwise stated. We define two subsets 
of sensors, Sa and Sf,, as follows. A sensor sj is in Sa if and only if there is no sensor Sj with 
i < j such that their covering intervals are strictly overlapped (e.g., see Fig. [6]). A sensor Si is 
in Sb if and only if there is no sensor Sj with i < j such that their covering intervals are strictly 
overlapped. Let the indices of sensors in Sa be ai, 02, • • • , a„j and the indices of sensors in Sh be 
61, 62, ... , 6„2, from left to right. We claim ni = n2. To see this, consider the interval graph G in 
which the covering interval of each sensor is a vertex and two vertices are connected by an edge if 
their corresponding intervals are strictly overlapped. Observe that in each connected component 
of G, there is exactly one interval whose corresponding sensor is in Sa and there is exactly one 
interval whose corresponding sensor is in Sb, and vice versa. Thus, ni = 77-2, which is the number 
of connected components of G. Let m = ni = n2 < n. Further, it is easy to see that the covering 
intervals of both Oj and bi must be in the same connected component of G and ai <hi. Indeed, Oj 
(resp., hi) is the leftmost (resp., rightmost) sensor in the subset of sensors whose covering intervals 
are in the same connected component of G (see Fig. [6]). Note that Oj = hi is possible. Hence, G 
has m connected components. 

For each 1 < z < m, let Gi denote the connected component containing the covering intervals 
of ai and hi; with a little abuse of notation, we also use Gi to denote the subset of sensors whose 
covering intervals are in the connected component Gi. Clearly, Gi = {sj | < j < hi} (e.g., see 
Fig. [6|). We also call Gi a group of sensors. The groups Gi, G2, . . . , Gm form a partition of S. The 
sensor Sq- (resp., Sh^) is the leftmost (resp., rightmost) sensor in Gi. 

Lemma 14 For any two sensors Si and sj with i < j, if Si ^ Sb or sj ^ Sa, then A3(i,j) 7^ A*. 

Proof: Assume Si Sb. In the following, we prove that A3(i,j) cannot be A* for any i < j. 

Suppose Si is in the group Gk- Then i < hk due to Si ^ Sb- Further, the covering intervals of 
Si and Sj+i must be strictly overlapped (otherwise, would be in Sb). Assume to the contrary 
Xsiijj) = A*, which implies that case (c) in Lemma [7| holds. Thus, in the corresponding OPT, Si 
is moved to the right by the distance X3{i,j) and all sensors Si,Si+i, . . . ,Sj must be in attached 
positions. It is easy to see that the sensor Sj+i must move to the right by the distance A3(i, j) + 2r — 
{xi+i — Xi). Since the covering intervals of Sj and Sj+i are strictly overlapped, 2r — (xj+i — Xi) > 0. 
Therefore, the moving distance of s.j-|_i must be larger than that of Sj. Since the moving distance 
of Si is X3{i,j) = A*, we have contradiction. Hence, A3(i,j) cannot be A*. 

Assume Sj Sa- Then by a symmetric argument, we can prove X3{i,j) 7^ A* for any i < j. □ 
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By Lemma if A* G A3, then it can only be in the set A3 = {A3(i,j) \ i < j,i £ Sb,j € Sa}, 
and IA3I = 0{m?). Thus, we seek to order the elements of A3 into 0{m) sorted lists and each 
list contains 0{m) elements. One might be tempting to use the following way. Clearly, for each 
l<A;<m — 1, A3 contains X^isbf., Sa^) for all h = k + l,k + 2, . . . ,m, and hence one may simply put 
them into a list. However, such a list is not necessarily sorted. Specifically, for any two indices hi 
and /i2 with k + 1 < hi < h2 < m, either X3isb^,Sa^^) < Xsisbk^Sa^^) or Xaisb^^, Sa^^) > X3{sb^,Sa^^) 
is possible. Our approach relies on additional observations. Below, for simplicity of notation, we 
use A3(6fc,a/j) to refer to A3(sb^,Sa^). We first examine the value of each X3{bk,ah) in A3. 

By definition, we have X3{bk,ak+i) = (2:0^+1 — Xb^ — 2r)/2, and this is equal to half the length 
of the interval between the right extension of Sb,, and the left extension of Sq^^^, which we call a 
gap. Note that this gap is when the two sensors Sb^ and Sq^.^^ are in attached positions. For 
each 1 < k < m — 1, define = Xa^^^ — Xb^ — 2r, which is the length of the corresponding gap. 
Hence, A3(6fc,Ofc+i) = gkf^- Further, for each 1 < A; < m, we define the width of the group Gk as 
the length of the union interval of the covering intervals of the sensors in Gk, and define Ik as the 
sum of the lengths of the covering intervals of the sensors in Gk minus the width of Gk, i.e., Ik is 
equal to 2r{bk — + 1) minus the width of Gk- We then have the following observation. 

Observation 3 For every k with \ < k < m — 2, we have A3(6fc,a/i) = {Ylt=k 9t ~ Ylt=k+i 
for each h with k + 2 < h < m. 

Proof: By definition, X3{bk,ah) = [{xa^—Xb^.—2r) — 2r{ah—bk — l)]/2. It is easy to see that the value 
Xah~Xbk is equal to X^tfc Qt plus the sum of the widths of all groups Gk+i, Gk+2, ■ ■ ■ , Gh~i, and 
the value 2r{ah — 6^ — 1) is equal to the sum of the lengths of the covering intervals of the sensors in 
the union of the groups Gk+i, Gk+2, ■ ■ ■ , Gh-i- According to definitions of If for /c + 1 < t < /i — 1, 
the observation follows. □ 
The following lemma will be useful later. 

Lemma 15 For four indices ki,k2, hi, andh2, suppose ma.x{ki,k2} < min{/ii, /i2}; then X3{bk^,ah^) 
-Xz{bk^,ah2) = Xz{bk2,ah^) - Xz{bk2,ah2), and consequently, X3{bk^,ah^) < Xz{bk^,ah2) if and only 
if X3ibk2,ah^) < Xz{bk2,ah2). 

Proof: Note that for every 1 < < m — 2, we have X3{bk,ah) = {^t=k 9t ~ Ylt=k+i^t)/'^ for 
k + 2 < h < m, and A3(6fc, ah) = gk/2 for h = k + 1. 

If hi = h2, then the lemma trivially follows since A3(6fcj,a/ij) = X3{bk^,ah^) and A3 (6^2)^/11) = 
Xj,{bk2, a/12). Thus we consider hi / /12, and only show the case with hi < /12 (the case with hi > /12 
is similar). By their definitions, we have X3{bk^,ah^) - X3{bk^,ah2) = {-J2t=h^ 9t + YJU'hl^t) l"^- 
Similarly, A3(6jt2, a/ij - A3(6fc2) «/i2) = (" Y!t^l 9t + Ylt=hl Hence, the lemma follows. □ 

Lemma [T5] implies that for any ki and k2 with 1 < ki < k2 < m—1, the sorted order of A3(6/t^, a^) 
for all t = ^2 + 1) ^2 + 2, . . . , m is the same as that of the list X3{bk2 , at) for t = ^2 + !> ^2 + 2, . . . , m 
in terms of the indices of at- This means that if we sort the values in the list A3(6i,0() for all 
t = 2,3, ... ,m, then for any 1 < k < m — 1, the sorted order of the list X3{bk,at) with all 
t = k + l,k + 2, . . . ,m is also obtained implicitly. Our "ordering" algorithm works as follows. 

We first explicitly compute the values X3{bi,at) for all t = 2, 3, . . . , m, which takes 0{m) time, 
and then sort them in 0(m log m) time. Let p be the permutation of 2,3, ... ,m such that the 
increasing sorted list of these X3{bi,at) values is A3(6i, ap(i)), A3(6i, ap(2)), . . . , X3{bi,ap^^_i-j). Note 
that the permutation p is immediately available once we obtain the above sorted list. For any 
1 < < m, we say the element X3{bk,ah) is valid ifk + l<h<m and is undefined otherwise. 



22 



By Lemma [T5l the valid elements in each list A3(6a:, cLp{i)), A3(6fc, ap(2))) • • • ) ^3ibk,(J'p{m-i)) are also 
sorted increasingly. Further, if we compute gi,g2, ■ ■ ■ , dm-i and li,l2, ■ ■ ■ Jm as well as their prefix 
sums in the preprocessing, then given the index of any valid element in the list, we can obtain its 
actual value in 0(1) time. Clearly, the preprocessing takes 0(n log n) time. Thus, we have ordered 
(implicitly) the elements of A3 into 0{m) sorted lists and each list has 0{m) elements. 

However, we are not done yet. Since eventually we will apply the binary search technique of 
Lemma [TU] to these sorted lists and the lists contain undefined elements, the algorithm may take 
an undefined element in such a list and use it in the decision procedure which is the algorithm for 
Lemma [8l But, the undefined elements do not have meaningful values. To resolve this, we assign 
(implicitly) to each undefined element an "appropriate" value, as follows. For each 1 < < m, let 
C{k) denote the list A3(6fc, ap(i)), X^i^k-, flp(2))) • • • ; -^3(6^, ap(m-i))- If 1 < A; < m, then the list C{k) 
has some undefined elements. For each undefined element, we (implicitly) assign an actual value 
to it such that the resulting new list is still sorted. The idea is inspired by Lemma [T5l We use the 
list C{1) as the reference list since all its elements are valid. Every other list C{k) has at least one 
valid element, for example, the element \z{bk-,o,k+i)- We compute explicitly the value A3(6fc,o/c+i) 
for each 1 < A; < m, in 0{m) time. For a list C{h) with 1 < k < m and any undefined element 
A3(^fc,«p(j)) in >C(fc), we assign to it (implicitly) the value A3(6fc,afc+i) + h{bi,ap^i)) - h{bi,ak+i) 
(note that all these three values have already been computed explicitly). The lemma below shows 
that the resulting new list jC{k) is still sorted increasingly with this value assignment scheme. 

Lemma 16 For any 1 < k < m, the list C{k) is still sorted increasingly after all its undefined 
elements are assigned values implicitly. 

Proof: Consider any k with 1 < A; < m, and any two indices i and j with 1 < i < j < m — 1. It is 
sufficient to prove \3ibk,ap^i)) < Xs{bk,ap(^j)). 

If both values are valid, then by Lemma [T5l the inequality holds. Otherwise, we assume 
A3(&fc, ap(j)) is undefined. After our value assignment, A3(6fc, ap(j)) = X3{bk,ak+i) + A3(6i, ap(j)) — 
A3(^i,afc+i)- Depending on whether A3(6i, Op(j)) is undefined, there are two cases. 

If A3 (61, ap(j)) is undefined, then we have A3 (5^, ap(j)) = A3 (6^, Ofc+i) + A3(6i, ap(j))- A3(&i, a^+i). 
Hence, A3(6fc, Op(j)) - A3(6fc, ap(j)) = A3(6i, ap(j)) - A3(&i, ap(i)) > due to j > i. If A3(6fc, ap(j)) is 
valid, then by LemmaUHl we have X3{bk,ap(^j)) = X3{bk,ak+i) + X3{bi,apQ)) - X3{bi,ak+i). Thus, 
X3ibk,ap(j)) - A3(6a:, ap(j)) = A3(6i, ap(j)) - A3(6i, ap(j)) > 0. 

Therefore, in both cases, we have A3(6fc,ap(j)) < A3 (61, ap(j)), which proves the lemma. □ 

In summary, in 0(n log n) time, we have (implicitly) ordered the elements of A3 into 0{m) sorted 
lists and each list has 0{m) elements such that every element in any list can be obtained in 0(1) 
time. Hence, Lemma [T3l is proved. We remark that assigning values to the undefined elements in 
A3 as above does no affect the correctness of our algorithm. Assigning values to undefined elements 
only makes our candidate set A for A* a little larger (by a constant factor), which obviously does 
not affect the algorithm correctness because the larger candidate set still contains A*. One might 
also see that the statement of Lemma [T3] (and thus Lemma [9|) is a little imprecise since we actually 
ordered only the elements in a subset A3 of A3 (not the entire set A3). 

3.3 The Special Uniform Case 

In this subsection, we consider the special uniform case in which all sensors are initially located on 
the barrier B = [0, L], i.e., < Xj < L for each 1 < i < n. We give an 0(n) time algorithm for it. 
Again, we assume A* > 0. 
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Clearly, Lemmas [6] and [7] still hold. Further, since all sensors are initially on B, in case (a) of 
Lemma [71 Si must be si. To see this, since si is initially located on B = [0, L], it is always the best 
to use si to cover the beginning portion of B due to the order preserving property. We omit the 
formal proof of this. Similarly, in case (b) of Lemma [71 sj must be Sn- We restate Lemma [7] below 
as a corollary for this special case. 

Corollary 2 If X* > 0, then in OPT, there exist a sequence of consecutive sensors Sj, Sj+i, . . . ,Sj 
with i < j such that they are in attached positions and one of the following three cases is true, (a) 
The sensor sj is moved to the left by the distance X* , i = 1, and yi = r. (h) The sensor Si is moved 
to the right by the distance X* , j = n, and yn = L — r . (c) For i ^ j (i.e., i < j), the sensor Si is 
moved to the right by the distance X* and the sensor Sj is moved to the left by the distance A*. 

For any 1 < i < j < n, we define Xs{i,j) in the same way as before, i.e., X3{i,j) = [xj—Xi—2r{j— 
i)]/2, which corresponds to case (c) of Corollary[2l For each I < j < n, define X'i{j) = xj +r — 2rj, 
which corresponds to case (a). Similarly, for each 1 < i < n, define X2{i) = L — 2r{n — i) — (xj + r), 
which corresponds to case (b). We still use A3 to denote the set of all X^{i,j) values. Define 
A'^ = {A'i(j) I 1 < J < n} and A'2 = {X'^ii) | 1 < i < n}. Let A' = A'^ U A'2 U A3. By Corollary [21 we 
have A* € A'. The following lemma is crucial to our algorithm. 

Lemma 17 The optimal value X* is the maximum value in A'. 

Proof: Let A' be the maximum value in A'. It suffices to show A* < A' and A' < A*. Since A* G A', 
A* < A' trivially holds. Below, we focus on proving A' < A*. 

Since A* > 0, A' > holds. Clearly, either A' G A'^^, or A' G A2, or A' G A3. Below we analyze 
these three cases. 

If A' G A'^, then suppose A' = A'^(j) for some j. Since A' > 0, we have < A' = A'^(j) = 
Xj + r — 2rj, and thus Xj — r > 2r{j — 1). Since all sensors are initially on the barrier B, xj < L 
holds. Hence, even if all sensors si, S2, . . . , Sj-i are somehow moved such that they are in attached 
positions to cover the sub-interval [0, 2r(j — 1)] of -B, the sub-interval [2r(j — 1) , Xj — r] of B is still 
not covered by any of the sensors si, S2, . . . , sj-i. By the order preserving property, to cover the 
sub-interval [2r{j — l),Xj — r], the best way is to move sj to the left such that the new position of 
Sj is at 2r(j — 1) + r (i.e., the sensors si, S2, ■ ■ ■ , Sj are in attached positions), for which the moving 
distance of Sj is exactly Xi{j). Therefore, the maximum sensor movement in any optimal solution 
has to be at least X[{j). Thus, A' = X[{j) < A*. 

If A' G A2, then the analysis is symmetric to the above case and we omit the details. 

When A' G A3, the analysis has a similar spirit and we briefly discuss it. Suppose A' = A3(i, j) = 
[xj — Xi — 2r{j — i)]/2 for some i < j. Since all sensors are initially on the barrier B, we have 
< Xi < Xj < L. Consider the sub-interval [xi + r, xj — r] oi B. Because A' > 0, we have 
Xj—Xi — 2r{j—i) > 0, and thus {xj—r) — {xi+r) > 2r{j —i — 1). This implies that even if we somehow 
move the sensors Sj+i, Si+2, . . . , Sj-i such that they are in attached positions inside [xi + r, xj — r], 
there are still points in [xi + r, xj — r] that are not covered by the sensors Sj+i, Sj+2; • • • , Sj-i. By 
the order preserving property, to cover the interval [xi + r, Xj — r], we have to use both Si and sj 
and the best way is to move Si to the right and move sj to the left by an equal distance so that 
all sensors Sj, Sj+i, . . . ,Sj are in attached positions, for which the moving distances of Si and sj are 
both X3{i,j) exactly. Therefore, the maximum sensor movement in any optimal solution has to be 
at least X3{i,j). Thus, A' = A3(i,j) < A*. 

In summary, in any case. A' < A* holds. The lemma thus follows. □ 
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Base on Lemma [T71 to compute A*, we only need to find the maximum value in A', which can 
be easily obtained in 0{n?) time by computing the set A' explicitly (note that |A'| = Q{n?)). Yet, 
we show below that we can find its maximum value in 0{n) time without computing A' explicitly. 

Lemma 18 The maximum value in A' can be computed in 0{n) time. 

Proof: Let Ai, A2, and A3 be the maximum values in the three sets A'^^, A2, and A3, respectively. 
It is sufficient to show how to compute Ai, A2, and A3 in 0{n) time. 

Both the sets A'^ and A2 can be computed explicitly in 0{n) time. Thus, we can find Ai and 
A2 in 0{n) time. Below, we focus on computing A3. 

Note that for each value A3(i, j) G A3 with i < j, we have X3{i,j) = [xj — Xi — 2r(j — i)]/2. For 
each 1 < t < n — 1, define zt = Xt+i —xt — 2r. Hence, X3{i,j) = {Ylt=i ^t)/2. This implies that find- 
ing the maximum value in A3 is equivalent to finding a consecutive subsequence of zi, Z2, . . . , Zn-i 
such that the sum of the subsequence is the maximum among all possible consecutive subsequences, 
which is an instance of the well studied maximum subsequence sum problem. This problem can be 
solved easily in 0{n) time. Specifically, we first compute all values zi, Z2, . . . , Zn-i, in 0{n) time. If 
all values are negative, then A3 is the maximum value divided by 2. Otherwise, we let Zq = 0, and for 
each 1 < t < n — 1, let = m.ax{zl_i ,0} + zt. It is not difficult to see that A3 = ^ • maxi<(<„_i{zj}. 
Hence, A3 can be computed in 0(n) time. The lemma thus follows. □ 

After A* is computed, we use the linear time decision algorithm for Lemma [8] to compute the 
destinations of all sensors such that the maximum sensor movement is at most A*. 

Theorem 5 The special uniform case of the BCLS problem is solvable in 0{n) time. 

4 The Simple Cycle Barrier Coverage 

In this section, we discuss the simple cycle barrier coverage problem and present an 0{n) time 
algorithm for it. Mehrandish [15j gave an O(n^) time algorithm by somehow generalizing the 
O(n^) time algorithm [8j for the uniform BCLS. 

In this problem, the target region R is on the plane enclosed by a simple cycle B that is the 
barrier we aim to cover. The sensors in S" = {si,S2, . . . ,Sn} are initially located on B and each 
sensor is allowed to move only on B (e.g., not allowed to move inside or outside R). All sensors in 

5 have the same range r. Here, the distance between any two points on B is not measured by their 
Euclidean distance in the plane but by their shortest distance along B. If a sensor is at a point p 
on B, then it covers all points of B whose distances to p are at most r. Suppose all sensors in S 
are initially ordered clockwise on B as specified by their indices. Our goal is to move the sensors 
along B to form a coverage of B such that the maximum sensor movement is minimized. 

Since S is a cycle here, a sensor is said to move clockwise or counterclockwise (instead of right 
or left). Let L be the length of B. Again, we assume L < 2nr (otherwise, it would not be possible 
to form a coverage of B). Since 5 is a cycle, if L < 2r, then every sensor by itself forms a coverage 
of B. Below, we assume L > 2r. Imagine that we pick a point po on the interval of B from Sn 
clockwise to si as the origin of B, and define the coordinate of each point p € B as the distance 
traversed as we move from pQ to p clockwise along B. Let the input coordinate of each sensor 
Si & S he Xi. Thus, we have < xi < X2 < • • • < x„ < L. Further, for each 1 < i < n, we let Si+n 
denote a duplication of the sensor Si with a coordinate = Xi + L, which actually refers to the 
position on B with the coordinate Xj. 
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Since all sensors have the same range, it is easy to see that there always exists an order preserving 
optimal solution OPT in which the sensors are ordered clockwise along B in the same order as that 
of their input indices. A formal proof for this is given in [15j . Again, let A* be the optimal moving 
distance. We can check whether A* = in 0{n) time. Below, we assume A* > 0. 

Actually, our algorithm considers a set of 2n sensors, S' = {si, S2, ■ ■ ■ , S2n}- Specifically, the 
algorithm determines a consecutive sequence of sensors, S'^j = {si, Sj+i, . . . ,Sj} C S' with 1 < z < 
j < i + n, and moves the sensors of S'^j to form a barrier coverage of B such that the maximum 
sensor movement is minimized. Clearly, for each sensor G S, at most one of Sk and its duplication 
Sk+n is in S'^j. In this simple cycle case, the definition of attached positions of the sensors is slightly 
different from that in the line segment case. In this case, we call the two endpoints of the covering 
interval of a sensor Sj its counterclockwise and clockwise extensions, such that when going clockwise 
from the counterclockwise extension to the clockwise extension, we move along the covering interval 
of Si. One sensor is always in attached position by itself. Two sensors are in attached positions if 
the clockwise extension of one sensor is at the same position as the counterclockwise extension of 
the other sensor. Note that unlike in the line segment case, if two sensors Si and sj are in attached 
positions, say, the clockwise extension of Si is at the same position as the counterclockwise extension 
of Sj, then since the sensors are on the cycle, it is possible that the clockwise extension of sj is 
in the interior of the covering interval of Sj (e.g., when L < 4r). Similarly, a sequence of sensors 
Si, Si+i, . . . ,Sj (with l<i<j<i + n) are in attached positions if the clockwise extension of st 
is at the same position as the counterclockwise extension of st+i for each i < t < j — 1 (and the 
clockwise extension of Sj may be in the interior of the covering interval of Sj). The next result is a 
corollary of Lemma [71 

Corollary 3 If X* > 0, then in OPT, there exist a sequence of sensors Si, Sj+i, . . . ,Sj in S' with 
l<:i<j<i + n such that they are in attached positions and the sensor Si is moved clockwise by 
the distance A* and the sensor sj is moved counterclockwise by the distance A*. 

For each pair of i and j with 1 < i < j < i + n, we define X{i,j) = [xj — Xi — 2r{j — i)]/2. Let A 
be the set of all such X{i,j) values. By Corollary El A* € A. The following result is similar to that 
of Lemma [T71 

Lemma 19 The optimal value A* is the maximum value in A. 

Proof: The proof is very similar to that for Lemma [T71 and we briefly discuss it below. 

Let A' be the maximum value in A. It is sufficient to show A* < A' and A' < A*. Due to A* E A, 
A* < A' trivially holds. Hence, we focus on proving A' < A*. Since A* > 0, we have A' > 0. 

Suppose A' = X{i,j) = [xj — Xi — 2r(j — i)]/2 for some i and j with l<i<j<i + n. Consider 
the clockwise interval [xj + r, xj — r] on B, i.e., the union of the points on B from Xi + r to xj — r 
clockwise. Since A' > 0, we have Xj — Xj — 2r(j — i)>0, and thus {xj — r) — (xj + r) > 2r(j — i — 1). 
This implies that even if we somehow move the sensors s^+i, Si+2; ■ ■ ■ , •Sj-i such that they are in 
attached positions inside [xi+r, Xj—r], there are still points in [xi + r, xj —r] that are not covered by 
the sensors Sj+i, Sj+2, . . . , By the order preserving property, to cover the interval [xi + r, Xj—r] 

on B, the best way is to move clockwise and move sj counterclockwise by the same distance, for 
which the moving distances of Si and sj are both X{i,j) exactly. Therefore, the maximum sensor 
movement in any optimal solution has to be at least X{i,j). Thus, A' = X{i,j) < A*. 

The lemma thus follows. □ 
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By using the same algorithm for Lemma [TSl we can find A* in A in 0(n) time. With the value 
A*, we can then easily compute an optimal solution (i.e., compute the destinations of all sensors) 
in 0{n) time, as follows. 

Suppose A* = A(i, j) G A for some i and j with l<i<j<i + n. In the case of i > n, we have 
j > n and let i = i — n and j = j — n. Thus, we still have A* = X{i,j) since X{i,j) = X{i — n,j — n) 
when i > n and j > n. Below, we assume 1 < i < n. Note that j > n is possible. 

First, we move Si clockwise by the distance A* and move sj counterclockwise by the same 
distance A*. Next, move all sensors Sj+i, Sj+2, . . . , Sj-i such that the sensors Sj, Sj+i, . . . ,Sj are in 
attached positions. Since A* is the maximum value in A by Lemma [T9l the above movements of 
the sensors Sj+i,Sj+2, . . . ,Sj-i are at most A*. Then, starting at the sensor Sj+i, we consider the 
other sensors Sj+i, Sj+2, • • • , Si-i of S clockwise along B, and move them to cover the portion of B 
that is not covered by the sensors Sj, Sj+i, . . . ,Sj. For this, we can view the remaining uncovered 
portion of i? as a line segment and apply the linear time greedy algorithm for Lemma [8] with the 
value A*. The overall running time is 0{n). 

Theorem 6 The simple cycle barrier coverage problem is solvable in 0{n) time. 

Note that the case L < 2nr is also discussed in |15) . where an issue of balance points appears. 
In the case L > 2nr that we consider, the issue does not exist because the entire cycle B must be 
covered by the sensors in the optimal solution. 
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